On the 14 C(a,7) reaction rate. 
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The 14 C(q,7) reaction rate at temperatures below 0.3 GK depends on the properties of two 
near threshold resonances in 18 0, the 1" at 6.198 MeV and the 3~ at 6.404 MeV. The a+ 14 C 
Asymptotic Normalization Coefficients (ANCs) for these resonances were determined using the a- 
transfer reactions 14 C( 7 Li,t) and 14 C( 6 Li,d) at sub-Coulomb energies. The 14 C(a,7) reaction rate 
at low temperatures has been evaluated. Implications of the new reaction rate on the evolution of 
accreting helium white dwarfs and on the nucleosynthesis of low mass stars during the asymptotic 
giant branch (AGB) phase are discussed. 
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PACS numbers: 21.10.Jx, 25. 70. Hi, 26.2Q.-f, 26.30.-k 



I. INTRODUCTION 

The 14 C(a,7) 18 reaction has been suggested to play 
an important role in several astrophysical environments. 
It was pointed out by Mitalas pj that electron capture 
on 14 N can produce 14 C in the degenerate helium core 
of low mass stars. The a-capture on 14 C then follows. 
It was suggested in that the 14 N(e~,^) 14 C(a,7) 18 
reaction (NCO reaction) may trigger the helium flash in 
the core of low mass stars earlier than the 3a reaction. 
The significance of the NCO reaction for the evolution 
of low mass stars was also considered by Spulak [3] who 
concluded that it is not effective compared to the 3a re- 
action. A contradicting result was obtained by Kaminisi 
and Arai Q who used the updated (but still very un- 
certain) 14 C(a,7) reaction rate and concluded that the 
NCO reaction is dominant for igniting the helium flash. 
While the influence of the NCO reaction on the evolution 
of low mass stars in the red giant phase is still uncertain 
it was firmly established by Nomoto and Sugimoto [f|, 
and later by Hashimoto, et al. that the NCO re- 
action may trigger the helium flash in accreting helium 
white dwarfs at a lower temperature and density than 
the 3a reaction. However, this conclusion is rather sen- 
sitive to the actual value of the 14 C(a,7) reaction rate at 
temperatures between 0.03 and 0.1 GK. 

The 14 C(a,7) reaction is also important for production 
of 19 F in asymptotic giant branch stars. Observations by 
Jorissen, et al., pf show enhanced fluorine abundance in 
the atmosphere of K, M, MS, S, SC, and C asymptotic gi- 
ant branch (AGB) stars. This finding indicates that 19 F 
is produced in these stars and its abundance can be used 
to constrain the properties of the AGB models. The path 
for 19 F production in an AGB star, as proposed by Joris- 
sen 0|, is rather complex, occurring in the He intershell 
of the AGB star. Neutrons from the s-process neutron 
generator reaction, 13 C(a,n), can be captured by 14 N, 
enhanced from the previous CNO cycle. This leads to 
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production of 14 C and protons through the 14 N(n,p) 14 C 
reaction. Then 18 is produced by the 14 C(a,7) reaction, 
or, alternatively by 14 N(a,7) 18 F(,5) 18 0. The presence of 
protons and 18 O isotopes in the He intershell simultane- 
ously allows for the 18 0(p,a) 15 N reaction to occur and 
subsequently 15 N(a,7) 19 F capture. Competing reactions 
that reduce the final abundance of 19 F are 19 F(a,p), 
18 0(a,7) 22 Ne and 15 N(p,a) 12 C. A detailed investigation 
of the 19 F production in AGB stars has been performed 
recently by Lugaro, et al., [1] and it was determined that 
the major uncertainties in the production of 19 F are as- 
sociated with the uncertainties in the 14 C (a, 7) 18 O and 
19 F(a,p) 22 Ne reaction rates. Addressing uncertainties of 
the 14 C(a,7) reaction rate at temperatures relevant for 
helium flashes in the accreting helium white dwarfs and 
nucleosynthesis in the asymptotic giant branch stars is 
the main goal of this work. 



Direct measurements of the 14 C(a,7) reaction cross 
section are only available for energies above 880 keV in 
cm. [§]. The relevant temperature range for accreting 
helium white dwarfs is between 0.03 and 0.1 GK and for 
19 F nucleosynthesis in AGB stars it is ~0.1 GK. The 
Gamow window for the 14 C(a,7) reaction at these tem- 
peratures corresponds to an energy range between 50 and 
250 keV. At these energies the cross section for the a cap- 
ture reaction on 14 C is too low to be measured directly 
and has to be extrapolated from the higher energy data. 
However, near threshold resonances in 18 O may have a 
very strong influence on this extrapolation. The level 
scheme of 18 O is shown in Fig. [TJ Three states with ex- 
citation energies close to the a decay threshold at 6.227 
MeV are known in 18 0. These are the 1" at 6.198 MeV, 
the 2- at 6.351 MeV and the 3" at 6.404 MeV. The 
2~ state at 6.351 MeV is an unnatural parity state and 
cannot contribute to the 14 C (0,7) reaction at any signif- 
icant level. The 1~ subthreshold state at 6.198 MeV can 
only contribute to the 14 C(a,7) reaction through its high 
energy tail, therefore its contribution is expected to be 
relevant only at the lowest energy. Finally, the 3~ state 
at 6.404 MeV is 177 keV above the a-decay threshold - 
right in the middle of the energy range of interest. 
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Figure 1: Partial level scheme of 18 0. 



The strength of a resonance is defined by 
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where T, T a and T 7 are the total, partial a and par- 
tial 7 widths and Jr, % and j p are spins of the reso- 
nance, the target and the projectile. At 177 keV the 
partial a-width is much smaller than the partial 7 width 
(r Q <C T 7 ), therefore the resonant a-capture cross sec- 
tion is determined entirely by the partial a-width with 
resonance strength calculated as ui^f w 7 x T a for the 3~ 
state at 6.404 MeV in 18 0. The partial a-width of the 
3~ state is not known. This leads to several orders of 
magnitude uncertainty of the 14 C(a,7) reaction rate at 
~ 0.1 GK. The subthreshold 1" state at 6.198 MeV also 
contributes to the reaction rate uncertainty at the lowest 
temperatures. 

Information on the partial a widths of resonances in 
O can be extracted using the direct a-transfer reac- 
tions 14 C( 6 Li,d) and 14 C( 7 Li,t). Usually these reactions 
are performed at energies of 5-15 MeV/ A and a spec- 
troscopic factors (S Q ) for the states of interest are de- 
termined from the Distorted Wave Born Approximation 
(DWBA) analysis of the reaction cross sections. The par- 
tial a width is then calculated as r a = S a x T SP , where 
Tsp is the a single particle width. 
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where Pi(kR) is a penetrability factor determined by the 
Coulomb regular and irregular functions Fi(k,R) and 
Gf,{k,R), /i is a reduced mass and k is a wave num- 
ber. R is a channel radius - R — ro(\ffii + v^Ai), where 
ro = 1.2 — 1.4 fm. The reduced width, j sp , in equation 
(j4]) corresponds to the reduced width in a square- well po- 
tential. The single particle width can be calculated more 
accurately using a realistic Woods-Saxon potential. As 
is well known, the spectroscopic factor extracted from 
an a-transfer reaction depends on the parameters of the 
optical potentials used to describe the wave functions of 
relative motion in the entrance and exit channels in the 
DWBA approach and on the shape of the form factor 
potentials and the number of nodes in the model core- 
cluster wave function. However, for astrophysical calcu- 
lations the asymptotic normalization coefficient (ANC) 
is needed instead of the spectroscopic factor, and if the 
transfer reaction is performed at sub-Coulomb energies, 
then the parametric dependence of the DWBA calcula- 
tions is greatly reduced. This approach has been applied 
before in Refs. Brune et al. Johnson et al. where 
the 12 C(a,7) and 13 C(a,n) reaction rates due to near 
threshold resonances in 16 and 17 were evaluated. In 
this work we measured the cross sections of the 14 C( 7 Li,t) 
and 14 C( 6 Li,d) a-transfer reactions at sub-Coulomb en- 
ergies and determined the ANCs of the 1" and 3~ states 
in 18 O. The discussion is structured as follows: the ex- 
perimental technique is discussed in the section |Ul the 
extraction of the ANCs and associated uncertainties are 
described in section \TU\ the new 14 C(a,7) reaction rate 
and it's astrophysical implications are discussed in sec- 
tion UVl and the conclusions are in section [V] 



II. EXPERIMENT 

The sub-Coulomb a-transfer reactions 6 Li( 14 C,d) and 
7 Li( 14 C,t) were performed at the John D. Fox Supercon- 
ducting Accelerator Laboratory at Florida State Univer- 
sity. The use of inverse kinematics, 14 C beam and 6 ' 7 Li 
target, has several advantages. First, it eliminates the 
background associated with ( 6,7 Li,d/t) reactions on the 
unavoidable 12 C admixture to the 14 C target. Second, 
inverse kinematics provides a boost to the light recoils, 
which is essential for the detection of the reaction prod- 
ucts, since the reaction is performed at low sub-Coulomb 
energies. 

Five silicon AE-E telescopes were used to detect 
the light recoils (deuterons/tritons) produced by the a- 
transfer reaction. The AE detectors ranged in thickness 
from 15 to 25 /zm, while the E detectors were 500-1000 
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Figure 2: (Color online) Particle identification spectra in AE- 
E silicon telescopes from 7 Li( 14 C,t) run. The x and y axes 
are energies in MeV deposited by the charged particles in the 
E and AE detectors, respectively. The intense group of pro- 
tons at 2 MeV (2.4 MeV total energy) corresponds to elastic 
scattering of 14 C by the hydrogen contamination in the 7 Li 
target. Tritons from the 7 Li( 14 C,t) reaction are highlighted. 



fim thick. Particle identification was performed using 
the standard AE-E technique. A sample of the 2D parti- 
cle identification plot is shown in Figure [21 where it can 
be seen that the separation of the different isotopes of 
hydrogen is sufficient for reliable identification. 

The lithium targets were prepared and transferred to 
the scattering chamber in a sealed container under vac- 
uum in order to prevent oxidation. The thickness of the 
lithium targets was ~ 20 iig/cm 2 and they were prepared 
on a Formvar backing. An accurate determination of 
the target thicknesses and detector solid angles was per- 
formed using elastic scattering of protons. The p+ 6,7 Li 
elastic scattering cross section at 95° with proton beam 
energy of 6.868 MeV was measured previously with a 3% 
accuracy in Bingham et al. [13], and this result allowed 
the product of the target thickness and the telescope solid 
angles to be determined by placing the telescopes at 95° 
with respect to the beam axis one by one and measuring 
p+ 6 ' 7 Li elastic scattering with each lithium target. 

For energies below the Coulomb barrier the angular 
distribution for a direct transfer reaction has a peak at 
back angles, 180° in the cm. frame which corresponds 
to a peak at 0° in the lab frame in inverse kinematics. 
Therefore, to measure the reaction cross section in the 
region of its maximum, the detectors were placed near 0° 
in the lab frame. The detector nearest to 0° was shielded 
from Rutherford scattering of the carbon beam on the 
lithium target with a 5 [im Havar foil in front of the 
telescope. This foil was thick enough to stop all of the 
carbon beam, but thin enough so that energy losses of 
the deuterons and tritons were small. 

The 14 C beam energy is determined by the condi- 
tion that the reaction cm. energy should be below the 



Figure 3: Spectrum of deuterons from the 6 Li( 14 C,d) reaction, 
measured by the 8° telescope at 8.8 MeV 14 C beam energy. 
The energy of the deuterons in the laboratory reference frame 
is shown on the x-axis. The 18 O states populated in this 
reaction are labeled. 



Coulomb barrier in both the entrance and exit channels. 
The 14 C beam energies of 8.8 MeV for the 6 Li target (2.64 
MeV in cm.) and 11.5 MeV for the 7 Li target (3.83 MeV 
in cm.) meet this requirement for population of near a 
threshold states in 18 O. 



III. ANALYSIS OF THE a-TRANSFER CROSS 
SECTIONS AND ANCS 

Spectra of deuterons from the 6 Li( 14 C,d) reaction and 
tritons from the 7 Li( 14 C,t) reaction at 8.8 MeV and 11.5 
MeV of 14 C beam energy, respectively, are shown in Figs. 
[3] and HI Peaks in the spectra correspond to the states 
or groups of states in 18 O and all of them can be iden- 
tified with the known states in 18 0. The near threshold 
states of astrophysical interest for this work are the V 
at 6.198 MeV and the 3" at 6.404 MeV. Peaks which 
correspond to these states appear at 2.5 and 3.0 MeV in 
Figs. [3]and[4l respectively. The experimental resolution 
of our setup in the cm. was 120 keV (FWHM), therefore 
these two states partially overlap. The following steps 
were performed in order to reduce the ambiguity in the 
extraction of cross sections from the measured spectra. 
The experimental resolution was accurately determined 
from the width of the peak which corresponds to the V 
state at 4.456 MeV. This state is separated from the other 
18 O states by more than 500 keV and is therefore well re- 
solved in our spectra. It was observed that the shape of 
the 1~ peak is slightly asymmetric with a characteristic 
"tail" at the lower energy end (see the inset in Figure |4]) 
probably due to radiation damage in the detectors. This 
detector response was then used in the further analysis 
by modeling it with a Gaussian function which has larger 
sigma below the centroid energy and lower sigma above 



4 




E(MeV) E(MeV) 



Figure 4: Spectrum of tritons from the 7 Li( 14 C,t) reaction, 
measured by the 8° telescope at 11.5 MeV 14 C beam energy. 
The energy of the tritons in the laboratory reference frame 
is shown on the x-axis. The 18 states populated in this 
reaction are labeled. The inset shows the double Gaussian fit 
of the well resolved 1~ state at 4.456 MeV which was used to 
determine the energy resolution of the experimental setup. 



the centroid energy (110 and 60 keV in the Lab. frame 
respectively) . All the peaks in the spectra are then fitted 
using this lineshape. 

Three states can contribute to the broad peak observed 
at ~ 3 MeV in the triton and ~ 2.5 MeV in the deuteron 
spectra. These are the aforementioned 1~ and 3~ at 
6.198 MeV and 6.404 MeV but also the 2" at 6.35 MeV. 
While this last state is an unnatural parity state and can- 
not be populated in a direct a-transfer reaction it can 
still contribute through the compound nucleus reaction 
mechanism. The magnitude of this contribution can be 
estimated using the 2~ state at 5.53 MeV as a benchmark. 
This state is clearly visible in the group of states at ~ 4.5 
MeV in the spectrum of tritons from the 7 Li( 14 C,t) reac- 
tion (see Figured]). Four states contribute to this peak. 
These are the 2" at 5.53 MeV, the 0+ at 5.336 MeV, the 
2+ at 5.255 MeV and the 3" at 5.098 MeV. The four 
state fit of the peak is shown in Figure [U There is also 
an unnatural parity 3 + state at 5.378 MeV, but its con- 
tribution is expected to be very small and since it has an 
excitation energy close to the + its possible contribution 
to the peak can be taken into account by overestimating 
the strength of the + state. Note that there are only 
four free parameters in this fit - the strengths of the four 
states. The positions of the states were fixed from the 
known values and the width is determined by the exper- 
imental shape as modeled above. In this way, the cross 
section for population of the 2~ state at 5.53 MeV was 
determined. Assuming that only the compound nucleus 
mechanism is responsible for population of the unnatural 
parity 2~ state at 5.53 MeV one can determine the cross 
section for population of the 2~ state at 6.35 MeV using 
the Hauser-Feshbach formalism and scaling it to the 2~ 



Figure 5: (Color online) Four state fit of the 4.5 MeV group of 
tritons from the 7 Li( 14 C,t) reaction. This fit was performed 
to determine the strength of the unnatural parity 2~ state at 
5.53 MeV which was used to evaluate the contribution of the 
compound nucleus mechanism to the measured cross sections. 



Table I: Optical-model parameters used in analysis of the 
14 C( 7 Li,t) and 14 C( 6 Li,d) reactions. 



V a R a W a WS R' a' R c Ref. 
Channel (MeV) (fm) (fm) MeV MeV (fm) (fm) (fm) 



14, 



1 C+ 7 Li 
18 0+t 

C+ 6 Li 250.0 3.26 0.65 



33.1 4.17 0.85 
130.0 3.43 0.72 



IS 



O+d 92.92 2.73 0.81 

v 3.53 0.6 

v 1.9 0.65 

t v 2.05 0.7 



14 C+a 
Q+d 

a- 



0.0 7.8 4.51 0.72 4.17 [13] 

8.0 0.0 4.19 0.80 3.43 [13] 

0.0 7.5 3.26 0.65 4.82 [14] 

0.0 2.5 3.65 0.71 3.4 [14] 

3.53 [15] 

1.9 [16] 

2.05 [16] 



a Form factor: Woods-Saxon. 

b Form factor: Woods-Saxon derivative, values in the table 
do not include a regular factor of 4: W=4W<j. 
v - varied to reproduce separation energies. 



state at 5.53 MeV. 

A sample of the three states fit to the peak which corre- 
sponds to ~ 6.3 MeV excitation in ls O is shown in Figure 
[6] The strengths of the 1~ and the 3~ states at 6.198 and 
6.404 MeV are the only two free parameters. Strength of 
the 2 _ state was fixed as described above. The same 
fitting procedure was performed for all telescopes. The 
angular distributions for the 6.198 MeV and 6.404 MeV 
states from the 6 Li( 14 C,d) and 7 Li( 14 C,t) reactions are 
shown in Figure [7] 

DWBA analysis of the angular distributions was per- 
formed using the code FRESCO (version FRES 2.4) [13]. 
The optical model potentials parameters for the DWBA 
calculations were adopted from [3] and Fortune and Ku- 
rath [H for the 6 Li( 14 C,d) and the 7 Li( 14 C,t) reactions 



respectively. The shape of the form factor potentials is 
from Kubo and Hirata [lB| and Buck et al. [HI]. These 
parameters are summarised in Table U It is important 
to note that both the X" and the 3~ states were treated 
as bound in the DWBA analysis, even though the 3 is 
actually unbound by 177 keV. This is a good approxi- 
mation since this state decays only by 7 emission. The 
dependence of the final result on the choice of the bind- 
ing energy will be discussed later in this section. The 
DWBA calculations were performed for beam energies 
at the center of the target, 11.42 MeV for the 7 Li target 
measurements and 8.7 MeV for the 6 Li target measure- 
ments. The ANCs for the 1~ and the 3~ states were 
extracted from the experimental data and the DWBA 
analysis. 

Following the results of [18] the ANC is defined in the 
following way. The radial overlap function of the bound 
state c consisting of two particles a + b (c=(ab)) can be 
approximated by a model wave function: 



(ab)lj 



(r) 



^(ab^j^nlj 



(r), 



(5) 



where ^nij( r ) ls tne bound state wave function for the 
relative motion of a and b, and S? ab ^ is the spectroscopic 
factor of the configuration (ab) with the corresponding 
quantum numbers I and j in state c. The asymptotic 
normalization coefficient C? b ^- defines the amplitude of 
the tail of the radial overlap function If^ij and at radii 
beyond nuclear interaction radius (r > Rpj) is given by 



L (ab)lj 



(r) - C, 



W- 



(ab)lj - 



r; c ,2 c + l 



/ 2 {2k ab r) 



(6) 



W_ ric j c +i/2 (2fc a 6T) is the Whittaker function describing 
the asymptotic behavior of the bound state wave function 
of two charged particles, k ab is the wave number of the 
bound state c (k a i, 
parameter r\ c - 



V2^a&£ c ) and r\ c is the Coulomb 
Z a Zb^ab/kab of the bound state c. The 
bound state wave function ^%xj ( r ) nas similar asymptotic 
behavior 



nlj ' 



2 (2k ab r) 
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where b c nl ^ is the single-particle ANC. From the equations 
©, I© and ©, 
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(8) 



Spectroscopic factors can be extracted directly from the 
experimental data 



da 

dQ ex p 



S1S2 



da 

dVlDWBA 



(9) 



where S\ is the a+ 14 C spectroscopic factor of the corre- 
sponding state in 18 O and S 2 is the a + d(t) spectroscopic 
factor of the ground state of 6 ' 7 ^Li. The squares of a+d(t) 
ANCs of the ground states of 6 Li and 7 Li are well known 
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Figure 6: (Color online) Fit of the triton spectrum from the 
7 Li( 14 C,t) reaction in the vicinity of 6.3 MeV of 18 O excita- 
tion energy. The black curve is the sum of the contributions 
from three states - the 1" at 6.198 MeV, the 2" at 6.35 MeV 
and the 3~ at 6.404 MeV. The red, blue and green lines rep- 
resent individual contributions of the 1~, 3~ and 2~ states 
respectively. The contribution of the 2~ state was not varied 
(see discussion in the text). 



and equal to 5.29 ±0.5 and 12.74 ± 1.1 fm _1 respectively 
Blokhintsev et al. Igamov and Yarmukhamedov plcf . 
Therefore, the ANCs for the 18 states populated in the 
( ' Li,d/t) reactions can be extracted. The R-matrix re- 
duced width of the state can then be determined from 
the ANC 1181 



7 C 2 



1 



Wl Vc> m/2 (2k ab R) 



2fi a b 



R 



b)lj> 



(10) 



where R is the channel radius. The partial a width 
r a of the state, which is the only missing parameter 
needed to determine the astrophysically relevant reso- 
nance strength (see discussion in section U]), is given by 
equations <[2j) and ([3]) in which 7^ is used instead of jg P . 

The DWBA cross section is proportional to the squares 
of the single particle ANCs determined by the form fac- 
tor potentials. Spectroscopic factors are inversely pro- 
portional to the single-particle ANCs (see eq. ([8])). This 
makes ANCs extracted from the experimental data using 
equations ([8]) and © insensitive to the details of the form 
factor potentials. Since the reaction is performed at sub- 
Coulomb energies the dependence of the DWBA cross 
section on the parameters of the optical model potentials 
is minimized. 

The partial a width of the 3" state at 6.404 MeV in 
18 O determined from the measured cross sections using 
the approach described above is T a = (7.8 ±2.7) x 10~ 14 
eV, and a detailed analysis of the parameter dependence 
of the result follows. The reduced width 7 Q is inversely 
proportional and the penetrability factor is directly pro- 
portional to the channel radius (see Eqs. (fTU|l and ((3])), so 
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that the explicit dependence of the partial a width on the 
channel radius cancels out. Only the implicit dependence 
through the kR argument of the Coulomb and Whittaker 
functions remains. The partial a width of the 3~ state 
calculated at the channel radius of 4 fm is 8.2 x 10" 14 eV 
and at the channel radius of 6 fm it is 7.6 x 10~ 14 eV. 
Therefore, the channel radius uncertainty contributes to 
the total error budget at the level of 8%. A channel ra- 
dius of 5.2 fm was adopted. 

The value of the spectroscopic factor extracted from 
the transfer reaction depends on the assumption of the 
number of radial nodes of the bound state wave function 
but the value of the corresponding ANC does not. The 
a+d(t) wave function which corresponds to the 6 ^ 7 ^Li 
ground state has one node (excluding the origin and in- 
finity). The number of nodes of the a+ C wave function 
for the 3~ state in 18 is not known. From the Talmi- 
Moshinski relation 2N + L = ^(2n.j +h), where the sum 
is over all nucleons in the cluster, the minimum number 
of nodes in the a cluster wave function of the 3~ state 
is 2. The DWBA calculations for 2, 3 and 4 node wave 
functions were performed producing r Q = 8.0, 7.9 and 
7.8 xlCP 14 eV respectively, so a 3% uncertainty is asso- 
ciated with the unknown number of nodes in the a cluster 
wave function of the 3" state at 6.404 MeV in 18 0. Three 
nodes were considered as the "best guess". 

The 3~ state is unbound by 177 keV but treated as 
a bound state with small binding energy. The choice of 
binding energy is arbitrary and may, in principle, influ- 
ence the final result. We performed calculations for sev- 
eral binding energies between 200 and 10 keV. While the 
dependence of the final result (value of r a ) on the specific 
value of the binding energy is very small there is a clear 
trend of decreasing T a with decreasing binding energy. 
Note, that the penetrability factor was always calculated 
at resonance energy, which is 177 keV above the decay 
threshold. The extrapolation of this trend into the pos- 
itive energy region suggests that the partial a width of 
a state at 177 keV in the cm. should be reduced by 
about 2% from the value obtained with 30 keV binding. 
The final value of T a — 7.8 x 10~ 14 eV takes this small 
correction into account. 

The 7 Li form factor potential was adopted from 
Kubo and Hirata [H] (Woods-Saxon with V=-91.2 MeV, 
R=2.05 fm and a=0.7fm). This potential has been shown 
to reproduce t+a scattering phase shifts at low energies 
as well as the a-particle binding energy in 7 Li (with one 
node and an i=l wave function). It was used also as the 
main coupling potential for the ( 7 Li,t) a transfer reac- 
tion. The dependence of T a on the specific choice of the 
Li form factor shape parameters was investigated. An 
increase of the radius parameter from 2.05 fm to 3.7 fm 
(sum of rms charge radius of a-particle and 3 H) increases 
r Q from 7.8 to 8.1xl0~ 14 eV. The potential depth is al- 
ways adjusted to reproduce the a+t binding energy with 
one node and an I = 1 wave function. Variation of the 
diffuseness parameter of this potential within reasonable 
limits (from 0.5 to 0.8 fm) has an even smaller influence 
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(a) 1" at 6.198 MeV from 6 Li( 14 C,d) 
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(b) 1" at 6.198 MeV from 7 Li( 14 C,t) 
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(c) 3" at 6.404 MeV from 7 Li( 14 C,t) 
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Figure 7: Angular distributions for the 1 at 6.198 MeV and 
3" at 6.404 MeV states from the 6 Li( 14 C,d) and 7 Li( 14 C,t) 
reactions measured at 8.8 and 11.5 MeV respectively. The 
solid lines correspond to the DWBA best fit. 



on the final value of T a (~ 2%). 

The dependence of the final result on the specific pa- 
rameters of the 14 C+a binding potential is also weak. 
The 18 O form factor potential was taken from Buck et al. 
[IU (Woods-Saxon with R=3.53 fm and a=0.6 fm). The 
potential depth was adjusted to produce the desired bind- 
ing energy. Variation of the radius parameter from 3.0 fm 
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to 4.2 fm results in ~ 4% variation of the partial a width 
and variation of the diffuseness parameter from 0.5 to 
0.8 fm produces ~ 10% variation. Note that the DWBA 
cross section depends strongly on the specific choice of 
the form factor potential parameters and varies by an 
order of magnitude with the changes described above. 
However, the value of the ANC is not sensitive to these 
variations. 

The influence of the optical model potential param- 
eters on the final result was minimized by performing 
the a-transfer reaction at a sub-Coulomb energy. Yet, 
since the 14 C beam energy was close to the Coulomb bar- 
rier uncertainty in the parameters of the optical model 
potentials gives largest contribution to the final uncer- 
tainty of the T Q value of the 3~ state. Assessment of this 
uncertainty is somewhat arbitrary but was determined 
here by performing a Monte Carlo variation of all the 
nuclear potential parameters. This process resulted in 
a Gaussian distribution of DWBA cross sections with a 
standard deviation of 30%. The radius and diffuseness 
changes to the imaginary part of the 14 C+ 7 Li optical po- 
tential contributed the most to the uncertainty. DWBA 
calculations were performed with several different sets of 
14 C+ 7 Li potentials taken from O, [3, [5l[ and verified 
that all of these (very different) potentials give r Q well 
within the 30% uncertainty limit determined using the 
Monte Carlo approach. 

The statistical uncertainty for the 7 Li( 14 C,t) 18 0(3~) 
reaction data is 7%. The combined systematic uncer- 
tainty in the absolute cross section normalization origi- 
nating from the determination of the product of the tar- 
get thickness times the solid angle is also 7%. Finally, 
there is a 9% uncertaint y in the known value of the a + t 
ANC, 12.74±1.1 fm -1 |20j]. Collecting all the compo- 
nents of the error budget yields the final relative uncer- 
tainty of the T a value of the 3~ state at 6.404 MeV in 
18 determined from the 7 Li( 14 C,t) a-transfer reaction 
at sub-Coulomb energy of 35%. 

An ANC can only be reliably extracted from the trans- 
fer reaction data if reaction is peripheral. Performing 
transfer reaction at sub-Coulomb energy is an important 
factor in achieving this condition. Limiting radius at 
which transfer reaction is considered in DWBA calcula- 
tions it was verified, that 95% of the transfer reaction 
cross section comes from radius beyond 7 fm. This is far 
outside of the sum of nuclear interaction radius of 7 Li 
and 14 C indicating that reaction is highly peripheral. 

The determination of the ANC of the 3~ state at 6.404 
MeV in 18 was also attempted from the 6 Li( 14 C,d) re- 
action data. Unfortunately, statistics were rather low 
and only an upper limit could be set on the cross section 
from the 6 Li( 14 C,d) data. Therefore, only an upper limit 
of ~ 10 -13 eV can be determined for the r a of this state 
which is consistent with the T a = (7.8 ± 2.7) x 10~ 14 eV 
determined from the 7 Li( 14 C,t) data. 

The ANC for the 1" state at 6.198 MeV and the cor- 
responding uncertainty was determined using the same 
steps as outlined above. The only difference is that this 
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Figure 8: (Color online) S-factor of the 14 C(a,7) reaction as 
a function of energy. The blue dashed line is the S-factor due 
to direct capture (adopted from Gorres et al. [§]), the red 
solid line is the S-factor due to the sub-threshold 1~ state at 
6.198 MeV in 18 O and the black solid line is the sum of the 
two without interference. The red dotted line corresponds to 
constructive and the black dash-dotted line to destructive in- 
terference of the 1~ state at 6.198 MeV with the direct capture 
amplitude. 



state is bound by 29 keV with respect to a-decay, there- 
fore no correction to the dependence of the ANC on the 
assumed binding energy was necessary. The square of the 
Coulomb modified ANC for this state was extracted from 
both the 7 Li( 14 C,t) and the 6 Li( 14 C,d) data sets. The re- 
sulting values are 2.6±0.9 fm -1 and 3.0±1.0 fm -1 . The 
two values extracted from the different reactions agree 
within error. As before, the uncertainty is dominated by 
the dependence on the optical model potential parame- 
ters. The final result for the squared Coulomb modified 
ANC of the 1" state at 6.198 MeV in 18 0, determined 
by combining values from both the ( 6 Li,d) and ( 7 Li,t) 
measurements, is 2.8±0.7 fm -1 . 



IV. THE 14 C(a,7) REACTION RATE 

With the T Q for the 3~ state at 6.404 MeV and ANC 
for the 1" at 6.198 MeV in 18 0, the 14 C(a,7) reaction 
rate can be reliably extrapolated down to very low tem- 
peratures. At temperatures below 1 GK four major con- 
tributors to the 14 C(a,7) reaction rate can be identified. 
These are the direct capture (DC), resonance capture 
through the 4+ and 3" states at 7.12 and 6.404 MeV and 
the sub-threshold resonance capture through the 1~ state 
at 6.198 MeV. Direct 14 C(a,7) capture was measured by 
Gorres et al. [§] in the energy range from 1.14 to 2.33 
MeV and was extrapolated to lower energies using the 
following parametrization of the astrophysical S-factor 

S(E) =2.18- 1.60E + Q.82E 2 . (11) 
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Figure 9: (Color online) The 14 C(a,7) reaction rate due to 
resonant and non-resonant capture. Resonant capture due to 
the 4+ state at 7.12 MeV and 3" state at 6.404 MeV are shown 
as the black dash-dotted and solid red curves, respectively. 
Direct capture is the red dotted curve, capture due to the 
1~ sub-threshold resonance at 6.198 MeV is the blue dashed 
curve. 

The resonance strength of the 4 + state was also deter- 
mined by direct measurements in 0| (ory = 0.46 ± 0.08 
eV). The resonance strength of the 3~ state determined 
through the ANC in this work is ivy = 7 x r Q = 
(5.5 ±1.9) x 10~ 13 eV. The resonance reaction rate was 
calculated as in Fowler et al. [22] 

N A < av > R = 1.54 x 10 n (/iT 9 )- 3 / 2 X LOJX 

x exp (=Hf8S t ) , (12) 

where /i is a reduced mass and T 9 is the temperature in 
GK. 

The astrophysical S-factor for the 1~ sub-threshold res- 
onance was calculated from the Coulomb modified ANC 
using the approach outlined in [18| as 

S{E) « (2L + 1)tt 2 4 l c T 7 x |CV | 2 , (13) 
\i \Hi -\- e c ) 

where |Cx- | 2 is the square of the Coulomb modified ANC 
(ANC divided by the corresponding Gamma function, 
T(L + 1 + rj c )) of the 1~ state, k c is the wave number of 
the V state k c = y/2fie c , £ c is the binding energy, T 7 is 
the partial 7 width of the 1~ state calculated from the 
known mean lifetime, T 7 = (1.8 ± 0.3) x 10~ 7 MeV. The 
combined relative uncertainty of the S(E) value is 30%. 
The S-factor due to direct capture is compared to the S- 
factor due to the 1~ sub-threshold resonance in Figure [8l 
It is important to note that at energies below 50 keV the 
contribution to the S-factor from the 1~ sub-threshold 
resonance is larger than from the direct capture. Based 
on the results of [§] direct capture is mainly due to the 



p-wave. Therefore, the direct capture amplitude may in- 
terfere with the 1~ sub-threshold resonance amplitude. 
The S-factors due to constructive and destructive inter- 
ference of these amplitudes are shown in Fig. [8] as red 
dotted and black dash-dotted curves, respectively. Obvi- 
ously, interference significantly amplifies the importance 
of the 1~ state. Unfortunately, the sign of this interfer- 
ence is not known and both cases will have to be con- 
sidered. The reaction rates due to direct capture and 
capture to the sub-threshold 1~ state at 6.198 MeV in 
18 O were calculated by numerical integration of: 

x exp(— 2nri) exp (— £f ) dE. 

The reaction rates due to resonance and non-resonance 
capture are shown in Figure [9] as a function of temper- 
ature. Three temperature regions can be identified. At 
T>0.3 GK the 14 C(a,7) reaction rate is dominated by 
the 4 + state at 7.12 MeV. The 4 + resonance strength 
was measured directly 0, US] and the uncertainty of the 
14 C(a,7) reaction rate in this region is 17%. In the tem- 
perature range between 0.03 and 0.3 GK the 3~ state 
at 6.404 MeV dominates. This temperature range is of 
particular interest for helium accreting white dwarfs and 
AGB stars. The 3~ resonance strength was determined in 
this work with an uncertainty of 35%, which determines 
the new uncertainty for the 14 C (0,7) reaction rate in the 
temperature range of 0.03-0.3 GK. Below 0.03 GK, direct 
capture and capture due to the sub-threshold 1~ state at 
6.198 MeV dominate. While direct capture was extrap- 
olated from direct measurements in Gorres et al. [§] and 
the S-factor due to the 1~ sub-threshold resonance was 
determined in this work with an uncertainty of 30% the 
reaction rate at this low energy region is still uncertain by 
two orders of magnitude due to the unknown interference 
between the direct and the 1~ sub-threshold resonance 
capture amplitudes (see Figure [8]). 

The new 14 C(a,7) reaction rate may have the most 
profound effect on the evolution of accreting helium white 
dwarfs. It was suggested by Hashimoto et al. [(| that un- 
der certain conditions the 14 N(e~ ,t/) 14 C(a,7) 18 (NCO) 
reaction dominates over the 3a reaction and triggers a he- 
lium flash at earlier times and lower density values than 
the 3a process would otherwise. However, this suggestion 
is sensitive to the actual 14 C(a,7) reaction rate. Various 
accretion rates have been considered in Hashimoto et al. 
[fj and it was concluded that the effect of the NCO re- 
action is larger for faster accretion. Following the evo- 
lutionary path of the central density and temperature of 
the helium accreting white dwarf with the accretion rate 
of 10~ 8 Mq yr -1 it was found (see Fig. 5 in Hashimoto 
et al. Q) that 14 C burning is ignited when the central 
temperature reaches 0.066 GK while the 3a reaction does 
not begin until 0.080 GK. This prediction was based on 
the hypothetical 14 C(a,7) reaction rate which turned out 
to be a factor of 30 higher than the rate determined from 
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Figure 10: (Color online) Comparison of the new 14 C(a,7) 
reaction rate (red solid curve) in the temperature range most 
relevant for the helium accreting white dwarfs to the rate used 
by Hashimoto, et al., 0| (black dash-dotted curve). 

our measurements. (Note that at these temperatures the 
reaction rate is dominated by the resonance capture due 
to the 3~ state at 6.404 MeV). Comparison of the reac- 
tion rate in the temperature range of interest for helium 
accreting white dwarfs used by Hashimoto and the one 
determined from our experimental data is shown in Fig- 
ure [TUl The new, much lower reaction rate would change 
the NCO ignition temperature to higher values. An esti- 
mate of the new ignition temperature based on equality 
of the corresponding 14 C(a,7) reaction rates gives a new 
value of about 0.075 GK which is very close, but still 
lower than the ignition temperature for the 3a reaction. 
Obviously, detailed calculations with the new 14 C(a,7) 
reaction rate are necessary to know whether the NCO 
reaction has an effect on the evolution of the helium ac- 
creting white dwarf. However, it is clear that the new, 
significantly reduced reaction rate casts serious doubt on 
the effectiveness of the NCO reaction chain as a trigger 
for helium flashes compared to the 3a reaction. 

It was found in Lugaro et al. [8|] that uncertainties in 
the 14 C(a,7) and the 19 F(a,p) reaction rates are the main 
contributing factors to the uncertainty of the production 
of 19 F in AGB stars. The authors of Lugaro et al. [H] 
cite 5 orders of magnitude uncertainty in the 14 C (a, 7) 
reaction rate due to the unknown spectroscopic factor 
of the 3~ state at 6.404 MeV. This uncertainty is elim- 
inated now. Surprisingly, the "recommended" 14 C (a,-f) 
reaction rate used in Lugaro et al. [8], which is based 
on the assumption that the spectroscopic factor of the 
3 state is 0.01 is almost exactly (within experimental 
uncertainty of this work) the same as the reaction rate 
determined from our measurements. This fortuitous co- 
incidence means that the yields of 19 F calculated in Lu- 
garo et al. @ with the "recommended" 14 C(a,7) reaction 
rate are accurate and no new calculations are necessary. 
Of course, the uncertainty due to the 19 F(a,p) reaction 



rate still remains and will have to be addressed. 

Finally, with the uncertainty for the 14 C(a,7) reac- 
tion rate reduced to just 35% in the temperature range 
of interest for the core helium flashes in low-mass stars 
the debate over possible importance of the NCO reaction 
chain as a trigger can hopefully be settled. 



V. CONCLUSION 

The 14 C(a,7) reaction rate was studied using an indi- 
rect approach. The asymptotic normalization coefficients 
(ANCs) for the 1" state at 6.198 MeV and 3" at 6.404 
MeV in 18 O were measured in the a-transfer reactions 
( 7 Li,t) and ( 6 Li,d) on 14 C. The measurements were per- 
formed at sub-Coulomb energies to minimize the depen- 
dence of the final result on the optical model potential 
parameters. The extraction of ANCs instead of spec- 
troscopic factors significantly reduced the dependence on 
the form-factor potential parameters and the number of 
nodes assumed in the a- 14 C wave function. The reso- 
nance strength of the 3~ state was determined from the 
ANC to be uj-f = (5.5±1.9) x 10~ 13 eV. It was found that 
the resonance capture due to the 3~ state dominates the 
reaction rate in the temperature range of 0.03<Tg<0.3 
which is the most relevant temperature range for the 19 F 
nucleosynthesis in AGB stars, evolution of helium ac- 
creting white dwarfs and core helium flashes of low-mass 
stars. At temperatures below 0.03 GK the 14 C(a,7) re- 
action rate is determined by the interplay between the 
direct capture and capture due to the sub-threshold 1~ 
resonance at 6.198 MeV. The S-factor due to the 1" sub- 
threshold state was determined in this work from the 
ANC and the direct capture S-factor was suggested in 
[§]. In spite of this knowledge, the reaction rate at the 
temperatures below 0.03 GK is still uncertain by two or- 
ders of magnitude due to the unknown interference (con- 
structive or destructive) between the two amplitudes. 

The new 14 C(a,7) reaction rate is a factor of 30 lower 
than the one used in Therefore, the importance of 
the NCO chain as a trigger for helium flashes in helium 
accreting white dwarfs suggested in [H] is reduced, if not 
eliminated all together. The "recommended" 14 C (a, 7) 
reaction rate used in [§] for 19 F nucleosynthesis calcula- 
tion in AGB stars fortuitously coincides with that deter- 
mined from our experimental data and the uncertainty 
of the 19 F production in AGB stars associated with the 
14 C(a,7) reaction rate is now eliminated. We hope that 
the new information on the 14 C(a,7) reaction rate will 
help to resolve the question of whether or not the NCO 
reaction chain can trigger helium flashes in cores of low- 
mass stars before the 3a reaction does. 

Authors are grateful to Peter Hoeflich and Akram 
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knowledge support of National Science Foundation under 
grant number PHY-456463. 



10 



[1] R. Mitalas, Astrophys. J. 187, 155 (1974). 

[2] K. Kaminisi, K. Arai, and K. Yoshinaga, Prog, of Theor. 

Phys. 53, 1855 (1975). 
[3] R. G. Spulak, Astrophys. J. 235, 565 (1980). 
[4] K. Kaminisi and K. Arai, Phys. Rep. Kumamoto Univ. 

6, 13 (1983). 

[5] K. Nomoto and D. Sugimoto, Pub. of the Astr. Soc. 

Japan 29, 765 (1977). 
[6] M. A. Hashimoto, K. Nomoto, K. Arai, and K. Kaminisi, 

Astrophys. J. 307, 687 (1986). 
[7] A. Jorissen, V. V. Smith, and D. L. Lambert, Astron. 

and Astrophys. 261, 164 (1992). 
[8] M. Lugaro, C. Ugalde, A. I. Karakas, J. Gorres, M. Wi- 

escher, J. C. Lattanzio, and R. C. Cannon, Astrophys. J. 

615, 934 (2004). 
[9] J. Gorres, S. Graff, M. Wiescher, R. E. Azuma, C. A. 

Barnes, and T. R. Wang, Nucl. Phys. A 548, 414 (1992). 
[10] C. R. Brune, W. H. Geist, R. W. Kavanagh, and K. D. 

Veal, Phys.Rev.Lett. 83, 4025 (1999). 
[11] E. D. Johnson, G. V. Rogachev, A. M. 

Mukhamedzhanov, L. T. Baby, S. Brown, W. T. 

Cluff, A. M. Crisp, E. Diffenderfer, V. Z. Goldberg, 

B. W. Green, et al., Phys.Rev.Lett. 97, 192701 (2006). 
[12] H. G. Bingham, A. R. Zander, K. W. Kemper, and N. R. 



Fletcher, Nucl. Phys. A 173, 265 (1970). 
[13] H. T. Fortune and D. Kurath, Phys. Rev. C 18, 236 
(1978). 

[14] A. Cunsolo, A. Foti, G. Imme, G. Pappalardo, and 

G. Raciti, Phys. Rev. C 24, 476 (1981). 
[15] B. Buck, C. B. Dover, and J. P. Vary, Phys. Rev. C 11, 

1803 (1975). 

[16] K. I. Kubo and M. Hirata, Nucl. Phys. A 187, 186 (1972). 

[17] I. J. Thompson, Comp. Phys. Rep. 7, 167 (1988). 

[18] A. M. Mukhamedzhanov and R. E. Tribble, Phys. Rev. 

C 59, 3418 (1999). 
[19] L. D. Blokhintsev, V. I. Kukulin, A. A. Sakharuk, D. A. 

Savin, and E. V. Kuznetsova, Phys. Rev. C 48, 2390 

(1993). 

[20] S. B. Igamov and R. Yarmukhamedov, Nucl. Phys. A 

781, 247 (2007). 
[21] J. E. Poling, E. Norbeck, and R. R. Carlson, Phys. Rev. 

C 13, 648 (1976). 
[22] W. A. Fowler, G. R. Chaughlan, and B. A. Zimmerman, 

Ann. Rev. Astron. Astrophys. 13, 69 (1975). 
[23] M. Gai, R. Reddy, D. A. Bromley, J. W. Olness, and 

E. K. Warburton, Phys. Rev. C 36, 1256 (1987). 



